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We present an all-magnetic scheme for the assembly and study of magnetic dipoles within designed 
confinement profiles that are activated on micro-patterned permalloy films through a precessing magnetic 
field. Independent control over the confinement and dipolar interactions is achieved by tuning the strength 
and orientation of the revolving field. The technique is demonstrated with superparamagnetic microspheres 
field-driven to assemble into closely packed lattice sheets, quasi- ID and other planar structures expandable 
into dipolar arrays that mirror the patterned surface motifs. 

The confinement and manipulation of charge carriers within intricate low- dimensional structures through 
designed electric potentials have led to revolutionary advances in semiconductor devices 1 . In contrast, the 
magnetic analog of creating designable energy landscapes for confining magnetic dipoles with concomitant 
tunability of the relevant forces solely through magnetic fields has however remained elusive. Such field-based 
approaches could offer powerful means to probe mesoscopic magnetic systems characterized by length scales 
ranging from tens of nanometers to micrometers, forces extending from femto- to nano-newtons and times scales 
above a microsecond. 

Schemes to confine interacting magnetic dipoles of micrometer to centimeter scale have been proposed as early 
as in 1878 2 . However, existing approaches 2 " 7 confront challenges such as hard- wall confinements that prohibit 
regulation of the size of the dipole cluster, restriction of magnetic dipolar interaction to either repulsive or 
attractive, and lack of control on introducing and removing a specific number of dipoles (smaller than millimeter 
size) within the interacting landscape. Magnetic traps have also been developed for the confinement of dipolar 
atoms 8 " 13 ; however, the current- carrying coils utilized offer limited configurations on the spatial design and size of 
the trap. 

Here we present a straightforward mechanism that overcomes these challenges by activating a time- orbiting 
potential on a micro-patterned permalloy thin film with a precessing magnetic field. Independent tunability of the 
soft-confinement force and dipolar interactions (reversibly between attraction and repulsion) as well as in-situ 
control on the number of confined dipoles are obtained solely through tuning the magnetic field. Moreover, the 
lithographically defined confinement scheme allows complex designer landscapes of desired length scales to be 
readily fabricated. Demonstrated by the remote, field-driven assembly of dipoles consisting of fluid-borne micro- 
spheres, this all-magnetic scheme serves as a prototype to stabilize various low- dimensional spatial constructs 
realizable over wide length scales. These features enable fundamental studies such as artificial atoms 14 " 18 , nuc- 
leation 719 " 24 , jamming 25 " 27 and frustration 28 " 30 as well as applications in biomedicine 31 " 39 , material assembly 40-45 , 
photonics 5,46 " 49 , magnetic logic 50 " 52 , chip based devices 53 " 57 and atom traps 8 " 13 . 



Results 

The architectures of the confining potential for magnetic dipoles are defined by permalloy (Ni 0 .8Fe 0 .2) thin-film 
patterns imprinted on a silicon substrate through photo or electronbeam lithography (see Methods section). 
Figure 1(a) schematically illustrates a random distribution of superparamagnetic microspheres (8 /urn in dia- 
meter, UMC4F/9560, Bangs Laboratories) in solution on a permalloy disk pattern (210 /urn in diameter, 40 nm 
thick) that have, in the absence of magnetic fields, been drawn to the surface by gravity. An externally applied field 
H ext = (H Xi H yy H z ) transforms the microspheres into dipoles while the permalloy disk becomes uniformly 
magnetized along the (H x , H y , 0) direction. The dipoles interact with each other as well as fields emanating from 
the disk 58 60 , arranging themselves into configurations that minimize the total energy (7, that has primary 
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Figure 1 | Formation of a time-orbiting potential on a permalloy disk. 

(a) Illustration of magnetic microspheres randomly distributed atop a 
permalloy disk within fluid environment on a silicon substrate in the 
absence of external magnetic fields, (b) In the presence of H ext = (60, 0, 60) 
(X 10 3 /47i) A/m, shown from top to bottom are the calculated energy 
landscape ( U disk ), simulation, and experimental snapshot of the system in 
equilibrium. Contour interval on the energy landscape is 2000 k B T y where 
k B is Boltzmann's constant and T = 300 K the temperature. Projections of 
the microspheres on the surface are drawn in circles in the simulation to 
show levitation of microsphere chains along H ext . Scale bar in the snapshot 
is 100 fim. (c) Same as in (b) but with H ext precessing about the z-axis at 
20 Hz with in-plane component Hi = 60 (X 10 3 /47i) A/m and angle of 
precession 6 = 45°, leading to a time-averaged energy landscape, L/jisk- 

contributions from the dipole-disk (U disk ) and dipole-dipole (U dip ) 
interaction energies (see Supplementary Information). As shown in 
Fig. 1(b), in the presence of a fixed, oblique external field of H ext = 
(60, 0, 60) (X10 3 /4tt A/m [ = (60, 0, 60) Oe]) oriented between the 
+x and +z directions, each dipole on the plane above the pattern 
experiences an asymmetric energy landscape, U disk , where a strong 
energy barrier and trench exist at the disk edges facing the — x and 
+x directions respectively. The typical force experienced by a micro- 
sphere (lying within a circular area of radius 80% that of the disk) is 
less than 1 pN, which corresponds to a time scale of greater than 
1.5 s for the dipole to move a distance equal to its own diameter. 
Therefore, precessing H ext about the z-axis at 20 Hz (time period 
0.05 s) results in the microsphere motion being insensitive to the fast 
varying component of the energy profile. Instead, the dipoles on the 
patterned disk respond to a time- averaged energy landscape, in 
the form of a bowl-shaped confinement profile surrounded by a 
circular trench as shown in Fig. 1(c). The precessing field also ren- 
ders rotational symmetry to the time-averaged dipole-dipole inter- 
action (Udip), facilitating formation of an ordered planar cluster 
rather than aggregated chains within the confining landscape. As 
discussed below, the time- averaged forms of these interactions are 
utilized for the confinement and assembly of the microsphere 
(dipole) cluster. 

We first obtain a quantitative understanding on the interplay 
between the two primary forces - the confinement force, 
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Figure 2 | Quantification of the two primary forces in the system: time- 
averaged confinement force (Fdi S k) and dipolar force (Fa ip ). In the 
presence of a precessing H ext as in Fig. 1(c), the in-plane (x) component of 
(a) Fdi S k and (b) Fd ip (exerted by another microsphere located at x = 0) are 
plotted as a function of x (in units of the disk radius i^disk and microsphere 
radius R). Fj isk and Fj ip are defined as the maximum in-plane confinement 
and dipolar forces respectively, (c) Contours of constant Fj isk and Fj ip are 
plotted in the parameter space spanned by H x and 6. 



Fdisk = — VUdi S k> that arise from the dipole-disk interaction and the 
dipolar force, Fdi P = — VU dip . In the presence of a precessing H ext 
oriented at, for example, 0 = 45° from the z-axis with an in-plane 
component Hi = 60 (X 10 3 /47i) A/m, the confinement force experi- 
enced by a microsphere atop the disk is < 1 pN directed towards the 
center as shown in Fig. 2(a). An interesting property is the sole 
dependence of the confinement force on H l5 since Fdisk oc Hi. 
Figure 2(b) shows that under the same precessing field, the corres- 
ponding dipolar force between two microspheres on the surface is 
—0.2 pN (repulsive) and weakens with their separation (x) as 1/x 4 . 
Moreover, since F dip ocH 2 (2 cot 2 0— l), the dipolar force can be 
tuned from being repulsive (6 < 0 C ) to attractive (0 > 0 C ), crossing 

over at 6 = 6 C = cot _1 ^l j V^—54.7 0 . Separate and independent 

control on the confinement and dipolar forces can thus be obtained, 
as illustrated in Fig. 2(c), by tuning H x and 0 along contours of 
constant Hi or constant H\ (l cot 2 9 — l) . Derivations for the above, 
including the expressions: 
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are available as Supplementary Information. 
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The ability to tune the dipolar interaction (with 0) while the soft- 
confinement force is fixed (with Hi) allows the confined dipoles to 
assemble into various structures. An important parameter character- 
izing the dipole structure is the average center- to -center microsphere 
separation, (Ar), determined by Delaunay triangulation 61 which 
serves as a measure of the cluster size. As shown in the plot of 
Fig. 3, (Ar) is varied over a factor of two from 4R to 2R (in contact) 
as 0 is increased from 35° to 70°. In the repulsive regime of dipolar 
interaction (0 < 0 C ), the dipoles form an expanded cluster within a 
center- attracting potential (insets A, B, C, F and G of Fig. 3), resem- 
bling the two-dimensional version of J. J. Thomson's atomic "plum 
pudding model" 62 and lending itself as an unstudied form of artificial 
atom 14-18 . In the attractive regime (0 > 0 C ), the dipoles form a sheet of 
closely packed hexagonal lattice centered on the disk (insets D, E, I 
and J of Fig. 3). Due to the initial dipole configuration, occasional 
jamming - inability to form a single lattice - could occur (inset I and 
F2 of Fig. 3) that is related to gelations and glass transitions 25 " 27 in 
bulk systems. The transition from a 2D plum pudding to an ordered 
lattice occurs at 0 = 0 C ~ 54.7° as borne out in the simulation where 
the lowest order dipolar interaction vanishes. A small deviation of the 
observed transitioning angle to 0 C ~ 60° in the experiment could 
arise from the weakly repulsive hydrodynamic forces due to rota- 
tional motion of the microspheres in a precessing field. The existence 
of multiple forms of competing interactions could cause frustra- 
tion 28 " 30 , leading to formation of exotic structures such as the loosely 
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Figure 3 | Tuning the structure of a 45 -microsphere cluster confined in 
the time-orbiting potential. The mean separation between microspheres, 
(Ar), determined by the bonding lengths given by Delaunay triangulation 
(lower left inset in the plot) is plotted as a function of 6 while H l = 60 
(X 10 3 /47i) A/m. Snapshots of the equilibrium microsphere configurations 
at corresponding values of 6 are shown below for both simulation (insets 
A-E) and experiment (insets F-J). Insets Fl and F2 are snapshots taken at 
12.6 s and 50.0 s after the instant switching of 6 from 35° to 90°, while inset 
HI is an enlargement from inset H. Lines are drawn between microspheres 
with less than 10 /mi center-to-center separation ( — 1.25 times the 
diameter) in insets Fl, F2 and HI. Scale bars are 50 /mi. 



connected chains observed at 6 = 57.5° (insets H and HI of Fig. 3). In 
this case, competition sets in between the shorter- ranged dipolar 
attraction and longer- ranged hydrodynamic repulsion that scale 
with inter-particle distance (x) as 1/x 4 and 1/x 3 respectively 6 . The 
soft- confining force acts as a catalyst in assisting the dipoles to over- 
come the energy barrier between them, thereby facilitating chain 
connections. In addition to equilibrium structures, transient effects 
are observed by abruptly changing the dipolar interaction (within 
1 /is). For example, reorienting the precessing H ext from the repuls- 
ive to attractive regime (9 = 35° ^ 90°) leads to nucleation of the 
cluster 719 " 24 - a first order phase transition that is often difficult to 
observe at the molecular scale and to reversibly produce in the 
experiments - as shown in inset Fl of Fig. 3. Reversible control on 
the cluster structure is demonstrated in Videos SI and S2 of 
Supplementary Information, where 6 is instantly switched from 
90° to 35° (Video SI) and from 35° to 90° (Video S2) while H x is 
held fixed at 60 (X 10 3 /4tt) A/m. 

Crucial to the study of cluster structure is the control on the 
number of objects. However, such control for magnetic dipoles has 
only been achieved through direct physical placement into and out of 
the confinement for millimeter- and larger sized objects 2 ' 5 ' 6 or 
through the random compartmentalization for much smaller 
objects 4 ' 7 . In contrast, by application of a sequence of fixed magnetic 
fields which disrupt the time- orbiting potential, or a rotating mag- 
netic field (e.g. about the x- orjy-axis) which "rolls" the microspheres 
along targeted direction, a desired number of dipoles can be loaded 
into or unloaded from the confinement zone in the presented 
scheme. Examples of loading dipoles from 50 to 60, 70, 80, 90 and 
100 followed by unloading are demonstrated in Videos S3 and S4 of 
Supplementary Information. 

Discussion 

The time-orbiting bowl-shaped profile can be readily extended to 
other shaped confining potentials with surrounding barriers and 
trenches defined by permalloy patterns. Shown in Fig. 4(a-d) are 
examples of triangle, square, long rectangle and octagonal ring geo- 
metries that maintain the same independent control on the confining 
and dipolar forces through 0 and Hi. These energy landscapes unveil 





\G = 90°l 



(7 



0=35° 











m 





Figure 4 | Time-orbiting potentials in the shapes of (a) triangle, (b) 
square, (c) long rectangle, and (d) octagonal ring. Simulations (Sim.) and 
experimental (Exp.) snapshots of the system in equilibrium are shown 
below and to the lower-right of the energy landscape in each panel under a 
precessing H ext with Hi = 100 (X 10 3 /4ti) A/m and 6 = 90° or 35°. Contour 
intervals on the energy landscapes are 2000 k B T. Scale bars below the 
experimental snapshots are 100 /mi. 
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previously unexplored opportunities to study field- driven assembly 
of magnetic dipoles beyond confinements with cylindrical symmetry 
or simple topology. For instance, the dipole cluster expands to reflect 
different trap motif and symmetries under repulsive dipolar inter- 
action; the long rectangle serves as a quasi ID magnetic channel; the 
discrete clusters decorating the octagonal ring is a result of the rare 
competition among inter-particle attractive force and the confine- 
ment force - two forces that are typically cooperative. The intricate 
balance between the interplaying forces and geometry of the soft- 
confining potential leads to the emergence of a complex and design- 
specific behavior of the dipoles. 

The advantages of the all-magnetic approach presented here 
include easy, non- contact manipulation of the dipoles by weak 
external magnetic fields, and absence of intricate wiring patterns 
for addressing specific electrodes as well as screening and heating 
effects that are generally present in charge-based approaches. We 
expect the presented scheme to have useful applications, for 
instances, in biomedical devices 31,32 : control of local concentration 
and optimal mixing can be achieved by the expansion and collapse of 
the dipole cluster; attracting the dipoles from different directions to 
the center of the disk serves as a simultaneous, multi-directional 
force probe on a biological entity attached to the center; dipole clus- 
ters of tunable spacing can also be fixed as filters to sort out objects of 
different sizes within a flow channel. From a materials stand 
point 40 " 42,46 ' 47 , through multiplexing an array of patterned traps, a 
two-dimensional field-tunable photonic grating becomes feasible: 
it offers a length scale relevant at near-visible wavelength and tunable 
over a factor of two. As its electric counterpart (semiconductor 
devices) in revolutionizing the electronic age, the all-magnetic 
scheme presented may also impact the fast-growing fields of mag- 
netic logic devices. The patterned time-orbiting potential could serve 
as a lithographically defined virtual circuit or wave guide for spins or 
spin waves in spintronics 50 and magnonics 51 , while the ability to 
create clusters of various sizes serves as a strategy to build islands 
of nanomagnets that are dipole-coupled to produce logic opera- 
tions 52 . We further envision the confinement scheme as a potential 
improvement over conventional traps for magnetic dipolar 
atoms 8 " 13 . Given other means to confine the atoms to two dimen- 
sions, a 10 /dm wide circular cavity (inverted disk) on the permalloy 
film for example, is capable of confining diamagnetic atoms carrying 
1 Bohr magneton with an effective spring constant of 10" 15 N/m 
under H x = 100 (X 10 3 /4tt) A/m and 6 = 45°. Such a pattern-based 
confinement scheme would allow one to explore the effect of trap 
geometry on the wave function of the condensate atoms. 

In conclusion, we have demonstrated a novel technique based on 
magnetic thin-film patterns and a precessing field to create time- 
orbiting potentials of desired architectures with independent control 
on the confinement and dipolar forces. As advantages over conven- 
tional confinement schemes, these lithographically defined soft- 
confining potentials can be remotely activated, micrometer-scaled 
dipoles loaded or unloaded in-situ and pattern-coded clusters of 
tunable structures and sizes achieved. With the ability to reversibly 
and rapidly tune the interplaying forces, structural phenomena are 
easily reproduced that enable studies in artificial atoms, nucleation, 
jamming and frustration. We also envision broad applications of this 
all-magnetic prototype in biomedical devices, assembly of functional 
materials such as field-tunable photonic gratings, magnetic logic 
devices, and, when appropriately scaled, trapping of cold quantum 
gases. 

Methods 

Fabrication of permalloy thin-film patterns. Permalloy thin-film patterns are 
imprinted on a silicon substrate by standard photolithography or electron-beam 
lithography: For the disk pattern, two layers of photo resist (LOR3B and SI 805, 
MicroChem) are applied to the Si substrate through spin-coating at 2500 rpm 
followed by baking at 180°C and 1 10°C on a hotplate for the respective layers; the disk 
patterns are then exposed through an optical microscope (BH-2, Olympus) and 



developed (in MF-319, Dow). For other patterns, two layers of e-beam resist 
(methylmethacrylate and polymethyl methacrylate, MicroChem) are applied to the Si 
substrate through spin coating at 4500 rpm followed by baking on a 180°C hotplate; 
pattern areas are exposed at 125 /iC/cm 2 using a scanning electron microscope 
(Helios Nanolab 600, FEI Company) and developed using 1 : 3 methyl isobutyl 
ketoneisopropyl alcohol (MicroChem). A 40 nm layer of Nio. 8 Fe 0 .2 (permalloy) is 
then deposited onto the resist-patterned surface by magnetron sputtering (ATC 
Orion, AJA International) followed by lifting off unwanted areas with acetone and 
60°C bath of rc-methyl-2-pyrrolidone. The resulting permalloy patterns are protected 
from the environment by another layer of Si0 2 (100 ~ 200 nm) through spin coating 
(Silicafilm, Emulsitone) at 4000 rpm and baking at 180°C for the disk patterns, or 
through magnetron sputtering for other patterns. The spin-coated Si0 2 surface is 
treated with NaOH (0.1 M) for 10 min prior to experiment. The dimensions of 
various patterns are (in fim): circular disk (210 in diameter), triangle (260 in side 
length), square (150 X 150), rectangle (75 X 300) and octagonal ring (200 X 250 with 
a 50 X 100 hole). 

Microsphere solution. Superparamagnetic microspheres (UMC4F/9560, Bangs 
Laboratories) with diameter of 8 fim are diluted in de-ionized water with 0.01 ~ 
0.05% Tween-20 (Sigma-Aldrich) for use on spin-coated Si0 2 surface or with 0.02 ~ 
0.05% Triton X-100 (Sigma-Aldrich) for use on magnetron- sputtered Si0 2 surface, 
preventing surface adhesion. The diluted microsphere solution is then placed on the 
permalloy-patterned substrate and contained within a polydimethylsiloxane ring 
capped by a cover glass to eliminate evaporation and fluid flow. 

Magnetic field generation. Computer program (Lab VIEW, National Instruments) 
controls the current sent by the power supplies (BOP 20-10ML, Kepco) into the 
electromagnets and solenoid, providing three-dimensional control of the external 
magnetic field H ext = (H x , H y , H z ) applied on the sample. To prevent net rotation and 
drifting of the microsphere cluster when a precessing field about the z-axis is present, 
the rotation of the in-plane field is programmed to cycle through the following 
sequence of angles at a speed of 20 revolutions per second: 180°, —360°, 180°, —360°, 
360°, - 180°, 360° and - 180°. Note that the time period is 0.3 s per cycle. 
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